Study of CP violation in D ->■ VV decay at BES-III 
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In this paper, we intend to study the problem of CP violation in D meson by D — ¥ VV decay 
mode in which the T violating triple-product correlation is examined. That would undoubtedly 
be another excellent probe of New Physics beyond Standard Model. For the neutral D, we focus 
on direct CP violation without considering D° — D° oscillation. Experimentally, by a full angular 
analysis one may obtain such CP violating signals, and particularly it is worth mentioning that the 
upcoming large D data samples at BES-III in Beijing will provide a great opportunity to perform 
it. 
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The topic of CP violation in the D-meson sector 
has been the subject of extensive studies involving both 
charged and neutral D meson decays these years [ll4l2|. 
In a recent reference (Hf , they exploited the angular and 
quantum correlation in the D a D° pairs produced through 
the decay of the ?/>(3770) resonance at BES-III to inves- 
tigate CP violation. They build CP violating observ- 
ables in e+e" -> ^(3770) -> D°D° -> {V x V 2 ){VzV A ) (V 
denotes vector meson) to isolate specific New Physics ef- 
fects in the charm sector [l3| ■ 

In practice, one can also probe the CP violation in D 
meson by one D decay without considering the quantum 
correlation, which is the motivation of this paper. Among 
the various kinds of D decay modes, D — > V1V2 and sub- 
sequently decaying to two pseudoscalars for each vector 
meson is a particularly interesting one in the perspective 
of the copious kinematics of final state interaction (FSI). 



As the same case in B meson [14|, [15| , the new type of T 
violating signal involving so-called triple-product (TP) 
will emerge by comparing a pair of CP conjugate pro- 
cesses, where TP is composed of the momentum of one 
vector meson and two polarizations. Assuming CPT in- 
variance, T violating TP asymmetry is equivalent to CP 
violating. We shall see in this letter such TP asymme- 
tries are related to the helicity components in the an- 
gular distribution of the D —> V1V2 process. Moreover, 
performing a full angular analysis is feasible and realistic 
in experiment. 

On the other hand, TP asymmetry (in the latter sec- 
tion we sometimes also refer to TP asymmetry as CP 
asymmetry assuming CPT invariance) is a sensitive sig- 
nal of New Physics. As is well known that Standard 
Model (SM) predictions for CP violation in the charm 
sector are very small, thus any significant such signals 
would be exciting. Currently, the BES-III experiment 
is collecting data at ip(3770) peak. In the short future, 
lots of events of D decay will be accumulated, which will 



provide a great opportunity to perform a full angular 
analysis to further achieve a valuable information for the 
question discussed here. 

Let us first consider the process D(p) — > 
Vi(k, ei)V~2(q, £2), where the two vectors Vi,V2 are 
characterized as their four-momenta and polarizations 
(fe,ei) and (5,62), respectively. We can write the most 
general invariant amplitude as a sum of three terms that 
we will call s,d,p [1J-117|, 
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where mi (m.2) is the mass of V± (V2), and the scalar 
coefficients a, b and c are generally complex and can re- 
ceive contributions from several amplitudes with differ- 
ent phases. Thus, one can parameterize the coefficients 
as fl4j 
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where aj, bj and Cj are the modulus of their corre- 
sponding complex quantities, and 5j denotes strong phase 
(also called unitary phase in Ref. [Hj]), and <fij is the 
weak phase which is the necessary condition for occur- 
ring of CP violation on the basis of Cabibbo-Kobayashi- 
Maskawa (CKM) mechanism [HELl in the SM. Squaring 
the matrix element of Eq. ([TJ , one obtains 
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Next, using CPT invariance, the matrix element for 
the antiparticle decay D(p) — > V\{k, ei)V 2 (q, e 2 ) can be 
written as 



X = 



with 



ae\ ■ e 2 + (p 4) 

m\m 2 

m\ra 2 



~ a = E 



iS s j„-i4>sj 



5 = J> 



(4) 
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Note that CP operator leaves strong phases invariant 
and only changes the sign of weak phase. From Eq. (fTJ) 
and Eq. (j4|), we will find that the p wave amplitude in 
M. changes the sign comparing with A4, which will in- 
duce an interesting property between \M\ 2 and \M\ 2 . 
To be clear, we would square the matrix element for the 
antiparticle decay in Eq. ((4]): 



\M\ 2 = \a\ 2 \e\ -e 2 \ 2 + -^K* ' ' 4)\ 2 + -^\e a ^ 5 e* la e* 20 k^ 5 \ 2 + 2?^-{e\ ■ e* 2 )(k ■ e* 2 )(q ■ ej) (6) 

_ 2 lMo£l {el . e . )6 a/M £ . ^ fc _ 2 M^) (fc . £ * )(g . e . )e a^ e J aC . p , . 

m\m 2 H mfm 2 H 

For Z?° — » V1V2 decay, one can define an asymmetry At with the definite sign for the triple product (k ■ e* x e|) 



as 



.14] 



JV(fc ■ e| x e* > 0) - N(k ■ e| x e* < 0) 
*4r = j- T , (") 

total 



where the subscript T implies triple products and N denotes the corresponding number of events. Eq. ([7]) above is 
actually 

= r(fc ■ e* x e* > 0) - £g ■ ej x e* < 0) 
T T(k ■ e* x e| > 0) + T(k ■ e\ x e* > ' 

Similarly, for £)° — > V\V 2 decay, At can also be constructed as the same way. In |.M| 2 , a triple-product correlation 
arises from interference terms involving the p amplitude, and will be present if Im(ac*) (or 7m(&c*)) is non-zero. 
After a simple calculation by inserting Eq. (|2|) and Eq. ([5]), we see that 



At cc Im(ac*) = ^ a,^ sin[(0 s . i - <\> vi ) + (S si - S pj )} (9) 

i,j 

Note that a non-zero triple correlation does not necessarily imply CP violation, since final state interactions (FSI) can 
fake it, namely the strong phase can also produce non-zero At (or At) even the weak phases are zero. Yet comparing 
a triple correlation for CP conjugate transitions allows to distinguish genuine CP violation from FSI effects. Thus 
we obtain, 

\{At + At) a ^[I m ( ac *) ~ Im(ac*)] = a;Cj sin((j) si - (p pj )cos(5 si - 5 pj ) , (10) 
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and 

\{At - At) cc i[/m(ac*) + Im(ac*)] = ^ cucj cos(0 sl - <j> p j) sin(£ s j - S pj ) . (11) 

So far, a non-zero value At + A7- wu l be undoubtedly a clean signal of CP non-conservation, because there must 
be at least one non-zero weak phase 4>. Here, we also note that if there is only one amplitude contributing to each 
partial wave, one can simultaneously determine the strong phase difference and weak phase difference from Eq. (|10[) 
and Eq. (ITTjl . 



Experimentally, one can perform a full angular analysis where 
to obtain the above TP asymmetry information because 
the complex coefficients a, b, c are related to the helicity 
amplitudes Ao, A\\, A± as discussed in Refs. flrl Eol]. 

Aq = —ax — b(x 2 — l), 
A\\ = V2a, 

Aj_ = y/2cy/x 2 - 1 , (12) 
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where wid is the mass of D meson. 

Now we turn to the full angular dependence of process 
D -> ViV 2 -> (PiP 2 ){P 3 P4 1 ) with P pseudoscalar, after 
some algebra one can get 



^ F n J± oc isin 2 6»i sin 2 6 2 cos 2 <f>\ Am | 2 + I sin 2 0i sin 2 6> 2 sin 2 </>|AjJ 2 + cos 2 6i cos 2 6> 2 |A | 

d cos Old cos 02dip 2 11 2 



sin 2 6i sin 2 2 sin 2(f>Im(A±A 



\/2 \/2 
H; — sin 2^i sin 26*2 cos 0i?e(A|| Aq) + — sin 20i sin 26> 2 sin 0Im(Aj_Ao) , (14) 
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where we have introduced A± with definite odd CP 
eigenvalue and the CP even partners Aq, Am via 
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(15) 



with Ax x \ 2 denoting the helicity mode of two vector 
mesons. 0[s (i=l,2) are the angles between the direc- 
tion of motion of one of the Vi,2 — > PP pseudoscalar 
final states and the inverse direction of motion of the 
D meson as measured in the Vi j2 rest frame, <j> is the 
angle between the two decay plane of vector mesons in 
the D rest frame. Figure [T] illustrates the decay kine- 
matics of the process D — > V{V 2 — > (P1P2XP3P4) in the 
rest frame of Vi ;2 . Eq. (|14p is consistent with that from 
Refs. |H[22|. 

As discussed previously, the TP asymmetry is 
connected with the angular dependence, combining 
Eqs. (JTDJl, CI]), GH), one can define the following T-odd 
quantities [14j, [l_5[ , 



A r ^ 



/m(A x A5) 



|A | 2 + |AJ2 + |A| 



(16) 




FIG. 1: Illustrative plot for the decay kinematics of process 
D -> V1V2 (P1P2XP3P4) in the rest frame of Vi, 2 . 



and 

Jm(Aj_A*|) 
Ar ~ |A | 2 + |A ± |2 + |A|||2' 

thus we will derive the CP violating observables, 
A = I(A r + A r ) 



(17) 



(18) 
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Im{A ± A%) 



Im(A ± A* ) 



+ |Aj 



A„ 



|A | 2 + |A. 



I4ll 



and 



A' = -(ArU + Ar 1 



(19) 



If- 

2V|A 



7ra(Aj_A 



\\> 



7m(A x Af 



2V|Ao| 2 + |Aj_p + |A|||2 lAoP + IA^P + IAnI 
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Before this study, there had been attempts to study 
CP violation of D meson via T violating TP corre- 
lation in theoretical viewpoint that differs from our 
method 12|, |23|, l24j . But the only reported experimen- 
tal search for T-odd asymmetries is from FOCUS in 
the D° -» K+K-ir+ir- and D+ -» K s K+tt + tt- de- 
cay modes [2Sj |. as listed in Table Q] No evidence for a T 
asymmetry is observed. The large BES-HI data sample 
is expected to provide enhanced sensitivity to possible T 
violating asymmetries. 





Decay mode 


A(%) 


D° 


-s- K+K--K+TV- 


1.0 ±5.7 ±3.7 


D+ 


->• KsK+TT + TT- 


2.3 ± 6.2 ± 2.2 


D + s 


-> KsK+it+iy- 


-3.6 ±6.7 ±2.3 



At last, we consider the potential sensitivity on the CP 
violating observables A and A' at BES-III. From Eq. (J7J , 
for a small asymmetry, there is a general result that 
it's error is approximately estimated as l/y/N tota i, where 
N total is the total number of events observed. At BES- 
III, with an integrated luminosity of 20 fb" 1 at ^(3770) 
peak, about 72 x 10 6 D°D° pairs will be collected with 
four year's running (2(| [2t]]. Table ITT1 lists some promis- 
ing channels to search for T asymmetry for both neutral 
and charged D decays and the corresponding expected 
statistical errors are estimated. The projected efficien- 
cies are extracted from Ref. 



obtained from Ref. 2S 



261 ] and branching ratios are 



TABLE I: T violating asymmetries in D meson decays from 
the FOCUS experiment 



VV Br (%) Eg. (e) Expected errors 



p°p°- 


* (7T + 7T )(n + TY ) 


0.18 


0.74 


0.004 


K*°p° - 


-¥ {K-ty+)(tt + ty-) 


1.08 


0.68 


0.002 


P°0-> 


(ir+ir-)(K + K-) 


0.14 


0.26 


0.006 


p+p- 


-> (7r+7r°)(7r-7r°) 


0.6* 


0.55 


0.002 


K*+K*~ 


->■ (K+-K )(K-n°) 


0.08* 


0.55 


0.006 


K*°K*° 


-J- (K + TT-)(K-TV+) 


0.048 


0.62 


0.002 


K*°p+ 


-)■ (K-ir + )(Tv+ir°) 


1.33 


0.59 


0.001 



TABLE II: The promising (VV) modes with large branching fractions, efficiencies and expected errors on the T asymmetry: 
the corresponding expected errors are estimated by assuming 20 fb" 1 data at ?/>(3770) peak at BES-III; the branching fractions 
with asterisk are estimated according to Refs. [29h32| . The last row is from D + decay. 



In tabic HH the branching fractions with asterisks have 
not been measured yet, but some estimates combining 
naive factorization and models for FSI are available from 
Refs. [HrHi]. Note that in table HI the estimated effi- 
ciencies are average value for the various partial waves by 
assuming that the magnitude of the longitudinal polar- 
ization is half of the decay rate. In the future, a careful 
measurements at BES-III about the efficiency and frac- 
tion for each partial wave are suggested. A more realistic 
analysis requires a likelihood fit to the full angular depen- 
dence of the D -> ViV^ -> (PiP 2 )(P3Pa) mode. System- 
atics will arise from the mis-reconstruction as V1V2 of the 
events that actually come from other resonances or non- 
resonance D — > P1P2P3P4 background contributions. In 
view of the sizable width of the vector resonances, we ex- 
pect that these systematics will dominate the final result. 
Their precise estimate in the BES-III experiment is be- 
yond the scope of this paper. However, as pointed out in 



Ref. [24| by Bigi, the four body decay of D -> P1P1P3P4 
both with and without intermediate states can all be used 
to probe T asymmetry in the frame work of TP. 

In conclusion, we studied the CP violation in D — > VV 
decay mode in which the T violating triple-product corre- 
lation is examined. That would undoubtedly be another 
excellent probe of new physics beyond the SM. The CP 
violating observables in connection with angular distribu- 
tion are constructed. For neutral D decays, we neglect 
the CP violation induced by D° — D° oscillation. Ex- 
perimentally, by doing a full angular analysis one may 
obtain such CP violating signals, and particularly it is 
worth mentioning that the upcoming large D data sample 
at BES-III will provide a great opportunity to perform 
it. The sensitivities for CP violating observables are esti- 
mated by assuming 20 _1 fb data-taking at -0(3770) peak 
at BES-III. 
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